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Abstract 

Energy is in general not conserved for mechanical nonholonomic systems with affine 
constraints. In this article we point out that, nevertheless, in certain cases, there is 
a modification of the energy that is conserved. Such a function coincides with the 
energy of the system relative to a different reference frame, in which the constraint 
is linear. After giving sufficient conditions for this to happen, we point out the 
role of symmetry in this mechanism. Lastly, we apply these ideas to prove that the 
motions of a heavy homogeneous solid sphere that rolls inside a convex surface of 
revolution in uniform rotation about its vertical figure axis, are (at least for certain 
parameter values and in open regions of the phase space) quasi-periodic on tori of 
dimension up to three. 
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1. Introduction 

This paper is ultimately addressed to the class of mechanical systems formed by a rigid 
body subject to the nonholonomic constraint of rolling without sliding on a surface which 
moves in a preassigned way. This type of nonholonomic constraints are affine, not linear, 
in the velocities. Consequently, even if the system is time-independent (which may easily 
be the case in presence of symmetries of body and surface, and if the latter moves at 
uniform speed) energy need not be conserved. In fact, while energy is conserved in time- 
independent nonholonomic systems with constraints that are linear in the velocities (see 
e.g. [22 122]), the same is not always true if the constraints are affine in the velocities (see 
e.g. [HE] and Section 3.5 below). A simple example of this situation is the classical 
system formed by a sphere that rolls without sliding on a table that rotates uniformly, 
studied by Pars [22, Neimark and Fufaev [22] and others, in which the energy is not 
conserved. 

*This work is part of the research projects Symmetries and integrability of nonholonomic mechanical 
systems of the University of Padova and PRIN Teorie geometriche e analitiche dei sistemi Hamiltoniani 
in dimensioni finite e infinite. 

^Universita di Padova, Dipartimento di Matematica, Via Trieste 63, 35121 Padova, Italy. Email: 
fassoOmath.unipd.it 

lUniversita di Padova, Dipartimento di Matematica, Via Trieste 63, 35121 Padova, Italy. Email: 
sanson@math.unipd.it 

§ Supported by the Research Project Symmetries and integrability of nonholonomic mechanical sys¬ 
tems of the University of Padova. 


1 



2 
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The starting point of this paper is the observation that for a nonholonomic system 
with affine constraints, even if the energy is not conserved, there might exist a modifi¬ 
cation of the energy—that may be interpreted as the energy of the system relative to a 
different reference frame and for this reason will be called moving energy —that is con¬ 
served. The reason underneath this fact is simply that, in a moving reference frame in 
which the surface is at rest, the constraint is linear. Therefore, if the system happens to 
be time-independent relatively to such a moving frame, its energy relative to that frame 
is conserved. And the condition of time-independence is easily verified in presence of 
symmetries. 

We study the existence of conserved moving energies in Section 3. For simplicity, 
instead of changing the reference frame with respect to which the system is described, 
we pass to a moving system of coordinates. After illustrating the mechanism on the well 
known example of a sphere on a turntable, we give sufficient conditions for the existence 
of a conserved moving energy (Theorem [TJ. Even though these conditions might appear 
very special, we relate them to the existence of symmetries: we assume that a group acts 
in configuration space and show that certain invariance properties of the system lead to 
a conserved moving energy (Theorem [2]). Interestingly, this conserved function is the 
sum of two functions that, at variance from what would happen in a holonomic system, 
are not conserved: the energy of the system and a component of the momentum map 
of the lifted action. There is here a connection with the failure of Noether theorem in 
nonholonomic mechanics, that we discuss. 

Lastly, in Section 4 we apply the results of Section 3 to the study of the system formed 
by a heavy homogeneous solid sphere that rolls without sliding inside a convex surface of 
revolution, which rotates around its (vertical) figure axis with constant speed fl. When 
= 0 this system is known to be integrable, with dynamics quasi-periodic on tori of 
dimension up to three [16] [27]. If fi ^ 0 the system is SO(3) x ^-invariant and, so far, it 
was only known that its four-dimensional reduced system admits two first integrals and 
a conserved measure, and thus that it is integrable by quadratures (6]. By exploiting 
the existence of an additional first integral given by a moving energy, much stronger 
integrability results can be obtained. Here we give a first, general result in this direction. 
Specifically, using essentially a continuity argument from the case fi = 0, we prove that, 
for small f2, there is an open nonempty subset of the reduced phase space in which 
the reduced dynamics is periodic, and correspondingly an open nonempty subset of the 
unreduced phase space in which the unreduced dynamics is quasi-periodic on tori of 
dimension up to three (Theorem [3]). 

Even though our primary interest is toward time-independent nonholonomic sys¬ 
tems, the need of considering time-dependent coordinates forces us to work in the time- 
dependent context. This somewhat complicates the notation. In the hope of keeping 
the complexity to a minimum, we adopt a Lagrangian description on the extended phase 
space of time-dependent mechanical nonholonomic systems, which is quickly described 
in Section 2. 

For general introductions to nonholonomic mechanics see e.g. EsunuziHniniisnBi; 
the time-dependent case is treated, using the formalism of jet bundles, in [20l l25l 1261 21] . 
Throughout the paper all manifolds and maps are assumed to be smooth, and all vector 
fields are assumed to be complete. 
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2. Nonholonomic systems with affine constraints 

2.1 Time-independent nonholonomic systems with affine constraints. First we 
briefly recall the time-independent case, mainly to fix some notation. The starting point 
is a Lagrangian system with n-dimensional configuration manifold Q and Lagrangian 
L : TQ —» R, that we assume to be regular; hence, in each set of local bundle coordinates 
(q,q) the matrix (q) is everywhere invertible. 

An affine nonholonomic constraint consists in the prescription that the kinematic 
states of the system belong to the fibers (Mo) q , q G Q, of an affine distribution Mo on 
Q, that we assume to have constant rank r > 1 and to be nonintegrable. Thus, there are 
a vector field £o on Q and a nonintegrable distribution T>o of constant rank r on Q such 
that 

(Afo )q = {v q G T q Q '. V q — £o (l?) € (®0 )q\ G Q ■ 

Of course, given Mo and X> 0 , Co is defined up to a section of Vq. V 0 will be said to 
be the distribution associated to Mo- The affine distribution Mo may also be regarded 
as a submanifold Mq of TQ, which is in fact an affine subbundle of TQ. We call Mq 
the constraint submanifold. The case of linear constraints is recovered when Co = 0; the 
constraint manifold is thus a linear subbundle of TQ. 

We assume that the nonholonomic constraint is ‘ideal’, that is, that it satisfies 
d’Alembert principle (see e.g. ! 2t EDI EE]) : at each q G Q, the set of reaction forces that 
the constraint can exert is (an appropriate jet extension of) the annihilator ( T>o)° q CT*Q 
of the fiber (T> o) q C T q Q of the distribution Vq associated to the constraint submanifold. 
It is well known that, under this assumption, there is a unique choice of the reaction 
force as a function 

Rl,m 0 '■ Mo —> Tq 

such that the restriction to Mq of Lagrange equations with the reaction force defines a 
vector field on M 0 (see e.g. mm • We denote this vector field on M 0 as Xl,q,m 0 and 
call it time-independent nonholonomic system , with affine constraints if Mq is an affine 
subbundle of TQ and with linear constraint if Mq is a linear subbundle of TQ. 

2.2 Definitions and notation for the time-dependent case. In order to consider 
time-dependent nonholonomic systems we pass to the extended phase space. In doing so, 
we need a number of definitions that we collect in this section. Let Q be an n-dimensional 
manifold. 

An m-dimensional extended submanifold M of TQ is an (m + l)-dimensional sub¬ 
manifold of the extended phase space TQ x R of the form 

M = {M t x{t} : t G R} . (1) 

Thus, for each t, M t is a submanifold of TQ of dimension m. The reason for the use of 
the term ‘extended’, instead of the perhaps more natural ‘time-dependent’, is that we 
need to treat both time-dependent and time-independent cases within the same context. 
We say that the extended submanifold M is time-independent if 

M = Mo x R 

for a given submanifold Mq of TQ, or equivalently if M t = Mq for all t, and that it is 
time-dependent otherwise. 

If all the Mt s in (jT|) are linear subbundles of TQ, then the extended submanifold M 
is an extended linear subbundle of TQ. If they are all affine subbundles, then M is an 
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extended affine subbundle of TQ. Obviously, we regard extended subbundles as special 
cases of extended affine subbundles. 

An extended distribution T> on Q is a distribution on Q x 1 with fibers 

= (Pt) q © { 0 } 

where, for each t, the (T>t) q are the fibers of a distribution T> t on Q. If all the distributions 
V t have rank r, then we say that T> has rank r. We say that V is nonintegrable if (some 
at least of) the distributions V t are nonintegrable. An extended distribution T> on Q of 
rank r generates an extended linear subbundle D = {D t x {f} : t £ R} of TQ of dimension 
n + r, with D t = {v q : q £ Q,v q £ (£> t ) g }, and vice versa. V is time-independent if 
T> t = V o for all t. 

An extended vector field on Q is a vector field £ on Q x R whose R-component is 
identically equal to 0, namely 


£(?>*) = 6 ( 9 ) + 0d t 

with each 6 a vector field on Q. £ is time-independent if 6 = £o for all t- 

If T> is an extended distribution on Q of rank r and £ is an extended vector field on Q, 
then 

M =V + £ 

is an affine extended distribution of rank r on Q. Thus M has fibers M.( q ,t) = (ffi > t)q+f,t(q), 
or M( q j) = (Mt) q © {0} with Mt = V t + 6- is called the extended distribution 
associated to M. M can be regarded in an obvious way as an extended affine subbundle 
M of TQ of dimension n + r. M is time-independent if so are £ and T>. 

Finally, a dynamical system on an extended submanifold M of TQ is a vector field 
on M whose R-component is identically equal to 1, that is 


X(v q ,t) = X t (v q ) + d t 


with each X t a vector field on M t . (The difference with respect an extended vector field 
is that now time does not stay constant, which is necessary for the dynamics). X is 
time-independent if X t = Xq for all t. 

2.3 Time-dependent nonholonomic systems. We start now from a Lagrangian sys¬ 
tem with ?r-dimensional configuration manifold Q and time-dependent regular Lagrangian 
L : TQ xR->R. The time dependency of the Lagrangian accounts, for instance, for the 
presence of time-dependent holonomic constraints))] 

We add now the nonholonomic constraint that its kinematical states belong to an 
extended affine subbundle M = {M t x {t} : t £ R} of TQ of dimension n + r, for 
some 1 < r < n, that we call the extended constraint submanifold. This extended affine 
subbundle corresponds to an extended affine distribution A4 = V + £ of rank r on Q. 
We assume that the associated extended distribution V is nonintegrable. 

The condition of ‘ideality’ of the constraint now means that, at each t and q , the 
set of reaction forces that the constraint can exert is (a jet extension of) the annihilator 

1 A large class of time-dependent holonomic constraints for systems of N material points can be mod¬ 

elled in this way. After the choice of a reference frame, that provides a (time-dependent) identification of 

physical 3-space with R 1 * 3 , a time-dependent holonomic constraint is given by a time-dependent embedding 
of a manifold Q into the configuration space (E 3 ) ;V of the unconstrained system, and the Lagrangian is 
the restriction of the Lagrangian of the unconstrained system to the resulting time-dependent, extended 
submanifold. 
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(2?f)° C T*Q of the fiber [V t ) q [201 ff) . and implies that there is a unique choice of the 
reaction force as a function 

Rl,m 

such that the restriction to M of Lagrange equations with the reaction force defines a 
dynamical system Xlqm on M. Here, T>° is the extended codistribution onQxl with 
fibers ( V t )° ® {0}. 

Definition. Let L : TQ x R —> R a regular Lagrangian and M an extended affine 
subbundle ofTQ. 

i. The dynamical system Xl^q^m on M is called the nonholonomic system with affine 
constraints (or, shortly, the nonholonomic system ) with Lagrangian L and extended 
constraint manifold M. 

ii. If M is an extended linear subbundle of TQ then we say that Xl,q,m has linear 
constraints. 

in. If L and M are time-independent, then we say that Xl^q^m is time-independent. 

In the time-independent case we will routinely identify Xl^q^m 0 and -Xl,q,Moxie > and, 
depending on the context, we will regard the Lagrangian L as defined either on TQ or 
on TQ x R. 

In (possibly time-dependent) bundle coordinates (q, q) in TQ, the fibers of the dis¬ 
tributions T>t on Q are the kernels of a (q, independent k x n matrix S(q,t) that has 
everywhere rank k, with k = n — r: 

C V t ) q = {q£T q Q : S(q,t)q = 0} . 

Thus, q £ (A it) q if and only if q = ft(q) + u with u £ ker S(q,t), that is, if and only if 
S(q, t)[q — £*(<?)] = 0. It follows that, for each t and q, the affine subspace {M-t) q of T q Q 
is described by 

S(q,t)q + s(q,t) = 0 (2) 

with s(q,t) = — S(q,t)£t{q ) G R fc . Of course, only ker S' is uniquely defined, not S, 
f and s. In coordinates, the annihilator of the fiber (T , t) g is the range of the matrix 
S(q,t) T , and the reaction force R.L,M(q,q,t ) £ range[S(g, t) T ]. 

2.4 Time-dependent diffeomorphisms and conjugation of nonholonomic sys¬ 
tems. In order to implement time-dependent coordinate changes, we need to consider 
(lifted) diffeomorphisms of the configuration space that depend on time, and use them 
to transform nonholonomic systems. (We may use now the expression ‘time-dependent’, 
instead of ‘extended’, because we will never need to consider ‘time-independent’ time- 
dependent change of coordinates and there will be no ambiguities). 

By a time-dependent diffeomorphism of a manifold U onto a manifold Q we mean a 
diffeomorphism C = ( Cq,C r ) : U x R —x Q x R. whose second component C K : U x R —> R 
is the identity between the R-factors. The first component Cq : U x M is a differentiable 
map, that in the sequel we denote Q. Thus, 

C = ( Q , id R ) : C/xR-MjxR, C(u, t) = ( Q(u , t ), t) 

and, for each t , the map 

Qt := Q(- ,t) : U^Q 

is a diffeomorphism (and “smoothly depends on f”). 
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With the identifications T(t/xl) ~ TU x TR and T(QxR) ~ TQ x TR, the tangent 
map TC : T{U x R) —x T{Q x R) can be seen as a diffeomorphism from TU x TR to 
TQx TR whose second component is the identity on the factor TR. Restricting TC to the 
unit tangent vector in the TR-factor gives the diffeomorphism TC |j =1 : TU x R x {1} —> 
TQ x R x {1}, that we regard as a diffeomorphism 

DC : TU xR-> TQ x R . 

Explicitly, if for all u £ U and f G R we write 

Qt{u) := G T Qt( „)Q, 

then 

DC(v u ,t) = ( T u Q t ■ v u + Qt(u), t ) G T Qt ( u) Q x R 

for all u G U, t G R and v u G T U U. Clearly, DC is a time-dependent diffeomorphism of 
TU onto TQ. 

In coordinates (u £ U, q £ Q) we will write Q' for ^ and Q for ^ and, for 
given t, Q' t = and Qt. = Thus, Q' and Q are defined on 17 x R 

while, for each t , Q' t and Q t are defined on U. With this notation, TC(u,u,t,t ) = 
( Q{u,t ), Q' t (u)u+ Q t (u)i , t , i) and 

DC(u,u,t ) = (Q t (u), Qj(u)u + Q t {u), t). (3) 

If C is a time-dependent diffeomorphism from U onto Q , then the pull back M := 
DC~ 1 (M) of an extended affine subbundle M = {M t x {f} : t £ R} of TQ is an extended 
affine subbundle of TU. In coordinates, if M is described by S(q,t)q + s(q,t) = 0 then 
M is described by 5(it, t)u + s(u, t) = 0 with 

S={SoC)Q', s = soC + (SoC)Q, (4) 

as is verified requiring that, for each t, (u, ii) £ M t if and only if (Q t (rt), Q' t (u)u+Q t (u )) G 
M t . 

The following fact is proven, in the time-independent case, in HU; the generalization 
to the time-dependent case is straightforward and we omit it. 

Proposition 1. Consider a nonholonomic system Xl^q^m cmd a time-dependent dif¬ 
feomorphism C from a manifold U onto Q. Then, the pull-back of under the 

restriction to M of DC coincides with the nonholonomic system Xp v ^ with L = Lo DC 
and M = DC^ 1 (M). 

3. Conservation of moving energy 

3.1 Example. We begin by illustrating the mechanism we have in mind on a sample 
system—the well known sphere on a turntable considered by Pars 23], Neimark and 
Fufaev |22t and several others, see e.g. urn- 

This system is formed by a homogeneous solid sphere constrained to roll without 
sliding on a table which, relatively to an inertial reference frame, rotates with constant 
rate f l around an axis orthogonal to it. In the mentioned references, and in all the other 
works we could find, the system is described with respect to the inertial reference frame, 


Francesco Fasso and Nicola Sansonetto 


7 


but is nevertheless time-independent. Let {0\x,y, z} be such a frame, and assume that 
the table lies in the sy-plane, and rotates about the z-axis. 

The configuration manifold Q is R 2 x SO(3) 9 (q,R) where q = (x,y) are the coor¬ 
dinates of the point of contact between sphere and table and R is the attitude matrix of 
the sphere. We identify TSO(3) and SO(3) x R 3 via right trivialization. The constraint 
manifold Mo is 8-dimensional and is diffeomorphic to R 2 x SO(3) x R 3 9 (q, R, ui), where 
ui = (u! x ,u) y ,u) z ) is the angular velocity in space of the sphere. Up to an inessential factor 
to, the mass of the sphere, the Lagrangian is 

L= l -\\q\\ 2 + l -ca 2 M\ 2 (5) 

where a is the radius of the sphere and ca 2 , with c > 0, its moment of inertia (divided 
by to). The condition of rolling without sliding is given by the affine constraint 

x = au)y — fly , y = —au) x + f lx . (6) 

The Lagrangian and the constraints are SO(3)-invariantH Reduction under this ac¬ 
tion consists merely in cutting away the SO(3) factor and produces a 5-dimensional 
reduced system on R 2 x l 3 3 (q, oj). The equations of motion of the reduced system are 
the two equations © and the three equations 

6j x = — (aoj v — fly) , Cj v = — (—aix x + fix), cJ z = 0 
a a 

where v = while the reconstruction equation is u! = RR T . It is elementary to show 
that the solutions of the reduced systems are periodic, with frequency u. 

In some of the quoted references, e.g. in [23], it is remarked that the reduced equations 
have the three independent first integrals 

ui z , aui x — vx , auiy — vy . 

However, the periodicity of a flow in a 5-dimensional phase space implies the existence 
of four, not just three, independent first integrals. To our knowledge (and surprise), this 
fact, and the existence of a fourth independent first integral, do not seem to have been 
noticed before. 

The obvious candidate for the missing first integral would seem to be the (projection 
to the reduced phase space of the) energy U/ J; m 0 of the nonholonomic system, which is 
the restriction to the constraint manifold Mq of the energy El of the Lagrangian L, that 
in this case coincides with L. However, energy is not conserved: the equations of motions 
give 

-^E l ,m 0 = acvfl(xuj y - y<jj x ). 

Nevertheless, a simple computation shows that the function 

E L ,M 0 -fl 2 (x 2 + y 2 ) + fLa(xuj x +yUy) (7) 

is a first integral of the nonholonomic system. Being SO(3)-invariant, this function is 
also a first integral of the reduced system, and it is independent of the previous three 
(except where x = y = 0). 


2 They are in fact invariant under an action of SO(3) X S 1 , but this is not used in the quoted references. 
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This additional first integral has a simple interpretation. In a system of time- 
dependent, rotating coordinates in which the table is at rest, the constraint of rolling 
without sliding on the table is linear, and the Lagrangian, which is the pull back L of L, 
is still time-independent. Therefore, the energy is now conserved and its push forward 
to the original coordinates is a first integral, that turns out to coincide with 0. 

The reason why the push forward of the energy in the rotating coordinates is different 
from the original energy (and may thus be conserved) is due to the time-dependency 
of the coordinate change. The Lagrangian in the rotating coordinates has the form 
L = L 2 + L\ + Lo, where the dependence of each Li on the velocities is homogeneous of 
degree The function L\ does not contribute to the energy E^ = L 2 — L 0 of L and 
the push forward of to the original coordinates differs from El by the push forward 
of L\. 

Remark. A completely similar situation is met in the system formed by a vertical disk 
constrained to roll without sliding on a uniformly rotating plane, considered in [131 . In 
that reference, the system is actually studied in a rotating frame, where the constraint is 
linear and the energy is conserved. However, the authors directly integrate the reduced 
equations of motion without noticing the conservation of energy. 

3.2 Moving energy and its conservation. First recall that the energy (or ‘Jacobi 
integral’) of a Lagrangian L : TQ x R —>• R is the function El : TQ xM->M given by 

E L {v q ,t) := ( p(v q ,t),Vq)q - L(v q ,t) Vq&Q, v q G T q Q , t€ R , 

where p is the momentum covector and ( , } q denotes the pairing between T*Q and T q Q. 
In coordinates, El = q • — L, where the dot denotes the scalar product in K”. 

Definition. Let Xl,q,m be a (either time-dependent or time-independent) nonholonomic 
system with affine constraints. 

i. The energy El,m ■ M —> R of Xl,q,m is the restriction of El to M: 

El,m '■= El\m ■ 

ii. If C : U xl —> Q xR is a time-dependent diffeomorphism, then the moving energy 
of Xl,q,m induced by C is the restriction Ef c M to M of the function 

e *l,c ■= EloDC o DC' 1 : TQ xl-)l. (8) 

Proposition 2. In the hypotheses of the above definition, 

E*L tC = E L -(p,QoC- 1 ). (9) 

Proof. The proof can be done in coordinates. If L = L o DC, then, from 0, L(u, u, t ) = 
L(Q t (u), Q! t {u)u + Qt{u), t). Thus Ei(u, u,t) = u ■ ||(u, u, t) - L{u , u, t ) = Q' t (u)u ■ 
§|( DC(u,u,t )) - L(DC(u,u,t )) = E L {DC(u,u,t )) - Q t (u) ■ ||( DC(u,u,t )). In 0 we 
have written Q 0 C -1 instead of QoDC -1 because Q is independent of the velocities. □ 

3 Clearly, L 2 may be interpreted as the kinetic energy, —Lq as the potential energy of the centrifugal 
force and —L\ as the generalized potential of the Coriolis force in a rotating, non-inertial reference 
frame in which the table is at rest. We prefer changing coordinates, instead of reference frames, since 
this exempts us from embedding the dependence on the choice of a reference frame into the theory, as 
e.g. in [20] . 
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The interest of considering a moving energy E £ c M resides in the fact that the func¬ 
tion Ef c differs from the energy El of the Lagrangian L by a term which is produced 
by the time-dependence of the diffeomorphism@ It is therefore possible that the function 
£/£ c M is a first integral even if El,m is not. We now formalize this possibility in the 
case of time-independent nonholonomic systems: 

Theorem 1. Consider a time-independent nonholonomic system with affine constraints 
Xl,q,m 0 x R an d a time-dependent diffeomorphism C from a manifold U to Q. Assume 
that: 

i. L o DC is independent oft. 

ii. Ef c is independent oft. 

in. DC^ 1 (Mq xl) is a time-independent extended linear subbundle of TU. 

Then, the moving energy Ef c MqxR is a time-independent first integral of Xl,q,m 0 xR- 

Proof. Hypothesis iii. means that DC~ 1 {Mq xR) = Mq x 1 for a fixed linear subbundle 
Mo of TU. By Proposition |I] DC conjugates Xl,q,m 0 xR to the nonholonomic system 
Xl u Mq xi w ith L = LoDC. By hypotheses i. and iii., X~ L fJ y f() xK has linear constraints 
and is time-independent. Therefore, the energy Ef ^ IoxR is a (time-independent) first 
integral of Xf v ^ R . It follows that its push-forward E^ ^ /oxR o (DC~ 1 \m 0 xr) is a first 
integral of Xl,q,m 0 x r- Since DC maps M 0 x R diffeomorphically onto M 0 x M, 

L,M 0 xR ° (DC |m 0 xr) = (^L,M 0 xR ° DC )|m 0 xR = E L ,C,M 0 xR ■ 

This proves that Ef c MoxR is a first integral of Xl,q,m 0 x r- Hypothesis ii. ensures that 
it is a time-independent function. □ 

We have stated Theorem 1 in terms of time-independent nonholonomic systems on 
the extended phase space so as to properly regard, in hypotheses i. and ii., the functions 
L and E’f c as defined onTQxR, even though constant on R. But if we identify functions 
on TQ and functions on TQ x R that are constant on R, then Theorem 1 states that, 
under hypotheses i., ii. and iii., a time-independent nonholonomic system with affine 
constrains Xl,q,m 0 has the time-independent first integral Ef c Mq (in fact, since Ef c 
is time-independent, Ef c MoxR and Ef c Mq may be identified). From now on, we will 
adopt this point of view. 

Remarks, (i) In time-dependent nonholonomic systems, either with linear or with affine 
constraints, the energy El,m is ordinarily time-dependent; even though it is not impos¬ 
sible that it is a (time-dependent) first integral, we do not consider this case because, in 
our opinion, from a dynamical point of view only time-independent first integrals are of 
interest. 

(ii) One might weaken the hypotheses of Theorem[l]in various ways, e.g. by requiring 
that only the restriction of Ef c to Mq x R be time-independent. However, the setting 
of Theorem [l] is sufficient for our application in Section [4] 

3.3 On the conditions of Theorem [1} The situation of Theorem [T] might appear very 
special. Our next goal, in Section 3.4, is to show that such a situation is instead easily 
verified in presence of symmetries. In order to gain some insight on this possibility, we 
begin by establishing conditions under which Ef c is time-independent and conditions 


4 A fact which is well known in the theory of time-dependent canonical transformations. 
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under which the extended constraint submanifold DC 1 (Mo x R) in the new coordinates 
is linear (even though possibly time-dependent). 

Proposition 3. Consider a time-independent nonholonomic system Xl,q,m 0 with affine 
constraints. Denote by £o + Vq the affine distribution on Q that corresponds to Mq, 
with £ 0 ® vector field on Q and T>q a distribution on Q. Consider a time-dependent 
diffeomorphism C = (Q, id K ) of Q to itself. 

i. Ef c is time-independent if and only if Q is the flow of a vector field Y on Q, and 
in that case 

E* Lfi = E L -(p,Y). 

ii. Assume that Q is the flow of a vector field Y on Q. Then, DC -1 (Mo x R) is an 
extended linear subbundle of TQ if and only if the vector field Y — £o is a section 
of Vo. 

Proof, (i) First note that Q o C~ 1 (q,t ) = Qt ° Qf (q) for all q and t. Let us write 
Y(q, t) = QtoQf 1 (q) and note that Y(q, t) € T q Q. Since L and El are time-independent, 
it follows from m that E* l c is time-independent if and only if (Y,p) is time-independent, 
that is, given that the momentum covector p does not depend on time, if and only if 

ft / BY \ 

gl( Y 'P) = ( apf 1 ) : re X R -> K 

vanishes. Since the Lagrangian L is regular, the map v q H > (p,v q ) q is a local diffeomor¬ 
phism for each q. Therefore, (^,p) vanishes identically in TQ x R if and only if ^ = 0. 
This shows that the time-independence of Ef c is equivalent to Qt o Qf 1 = Y with Y 
independent of t. But then Y is a vector field on Q and, since Q = Y o Q, Q is the flow 
of Y. 

(ii) The proof can be done in coordinates. Let M 0 be given by S(q)q + s(q) = 0. Then 
DC~ 1 (M 0 x R) is an affine subbundle of TQ that is described by S(u,t)u + s(u,t) = 0 
with S and s as in ©• Its linearity is equivalent to the vanishing ofs = soC + (S'oC)Q, 
that is, given that S and s are time-independent, s o C = s o Q and S' o C = S o Q, to the 
vanishing of s + S(Q t o Qf 1 ) = S(Y — £q)- □ 

This proposition suggests that, in order to obtain a time-independent conserved mov¬ 
ing energy, the time-dependent diffeomorphism C should be constructed as the flow of a 
vector field on Q that differs from the vector field £o by a section of the distribution Vq. 
The freedom in the choice of this section might then be used to try to make L o DC and 
DC~ 1 (Mq x R) time-independent. In the next section we will show that this is always 
possible if the system admits a symmetry group, with suitable properties, by choosing Y 
as an infinitesimal generator of the group action, that is, by choosing Q as the action of 
a one-parameter subgroup. 

3.4 Symmetry and conservation of moving energy. We consider now a time- 
independent nonholonomic system Xl,q,m 0 with affine constraints whose Lagrangian 
and constraint distribution are invariant- in a sense made precise in Hypotheses (HI) 
and (H2) below- under an action T : G x Q —> Q of a Lie group G on Q. As in 
Proposition [3l we denote by Afo = fo + Vo the affine distribution on Q corresponding 
to M 0 . 
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For each q G Q, we write as usual for g , q). We denote by : G x TQ —»• 

TQ the tangent lift of the action if, which is the action of G on TQ given by 

V* Q (v q ) = T q \\> g -v q 

(in coordinates, (q), A>' g (q)q) with ^' g = We make the following 

two hypotheses: 

(HI) L is invariant under namely 

L o Al T g Q = L VjeG 

(in coordinates, L($ g (q),V' g (q)q) = L(q,q) Mg,q,q). 

(H2) The distribution T>q is invariant under T, in the sense that 

(A))^ 9 (g) = Tq^g ■ (Do) q Mg £ G ,q G Q 

(we need not make any hypothesis on the nonhomogeneous term £o and on the 
invariance of Mq under the group action). 

Under these hypotheses, it is rather natural to try to build the time-dependent diffeo- 
morphism C = (Q,id R ) that leads to a conserved moving energy by choosing Q as a 
one-parameter subgroup of the action 'F. 

For 77 € 0 , the Lie algebra of G , denote by 

d 

Yr) := exp(tri)\t=0 

the infinitesimal generator of the action of the one-parameter subgroup generated by rj 
and by 

Jr) := {p, Yrj) 

the momentum map of the lifted action of the same one-parameter subgroup. The moving 
energy of Xl^q^m 0 relative of the time-dependent diffeomorphism C v = ( c f> y ’ I ,id R ), where 
: Q x K. — > Q is the flow of Y. q , is thus the restriction to Mq of the function 

Kc v =El-J v . ( 10 ) 

Theorem 2. Consider a time-independent nonholonomic system Xl^q^Mq with affine 
constraints and an action of a Lie group G on Q. Assume (HI), (H2) and 

(H3) rj € g is such that Y v — £0 is a section ofV 0 . 

Then, the moving energy Ef c Mq is a time-independent first integral of Xl^q^m 0 - 

Proof. Let Q be the flow of Y n . The conclusion follows if we show that the three hypothe¬ 
ses of Theorem [Q are satisfied with C = (Q, id R ). We may check them in coordinates. 

Hypothesis i. With this choice of Q , Q t = ’Fexp (tq) f° r all t- Hence, by hypothesis 
(HI), L(Q t (q), Q't(q)q) = L(q , q) for all q, q, t. Define L = L o DC. Since a vector field is 
invariant under its own flow, Q t = Y v o Q t = Q' t Y v . Thus 

L(q,q,t) = L(Q t (q), Q’ t (q)q + Qt(q)) = L(Q t (q), Q' t {q)[q + Y v (q)}) = L(q,q + Y v (q)) 

that shows that L is time-independent. 
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Hypothesis ii. Under hypothesis (H3), by item ii. of Proposition [3] M = DC~ 1 (M 0 x 
R) is an extended linear subbundle of TQ. We prove that this subbundle is time- 
independent. If Mo is described by S(q)q + s(q) = 0, then its associated distribution 
V 0 has fibers kerS'(g) and hypothesis (H2) is ker 5 , ('h s (g)) = ’I '' (q ) ker S(q) for all g and 
g, or 

ker [S'(Qt(g))] = Q' t (q) ker[5(g)] Vt,g. (11) 

In turn, the extended distribution M associated to M has fibers {Mt)q = ker[5(g, t)] 
with 

S(q,t) = S(Q t (q))Q' t (q), 

see ©■ Hence, using m and the fact that, if B is an invertible n x n matrix and S is 
a k x n matrix, then ker(5'l?) = B~ l ker S', we have 

ker S(q, t) = ker [S{Q t (q))Q' t (qj\ = Q((g) _1 ker S(Q t (g)) = kerS(g). 

Thus A4 t = T> o for all t and M is time-independent. 

Hypothesis in. This follows from item i. of Proposition [3] □ 

Remark. The condition that Y v — £o is a section of Vq implies that the orbits of the group 
action >I' T( 2 must be transversal to the constraint manifold. 

Example. The sphere on the turntable of Section 3.1 is an instance of the situation of 
Theorem 2. As in that section, we identify the tangent spaces to Q with R 5 via right- 
trivialization of TSO(3). The natural symmetry group of the problem is S 1 x SO(3) (that 
acts as in (fl^l) below), but for the sake of applying Theorem 2 we may consider only its 
subgroup G = S 1 x {!}, that acts on Q = R 2 x SO (3) as 

* e (q,R) = (H 9 q,H 9 R). 

Here H 9 is the matrix of the anticlockwise rotation of angle 6 around the 2 -axis and, 
with a small abuse, we identify vectors (x,y) € R 2 with vectors (x,y, 0) G R 3 . The 
(right-trivialized) infinitesimal generator of the action that corresponds to the Lie algebra 
element r/ G R is Y v = (— gy, px, 0, 0, q), and the corresponding momentum is J. v = 
rj(xy — yx + ca 2 to z ). The (right-trivialized) tangent lift of this action is 

T^ Q (g, R , q. u) = ( H e q , H e R , H g q, H 9 lo) 

and leaves the Lagrangian <[5|) invariant, as in (HI). The distribution Vq associated to 
the affine constraint © is given by x = aoj y , y = —aia x (that is, q = auj x e z with 
x the vector product in R 3 ) and is invariant under the action Te, as in (H2). Finally, 
the nonhomogeneous part of the constraint © is the vector field £ 0 = (—Uy, f lx, 0, 0,0) 
and Yq — £o = (0,0,0,0,12) lies in £> 0 , as in (H3). By Theorem 2, the moving energy 
El Cn m 0 = El,m 0 ~ Jq.\m 0 is conserved. Using © one verifies that this moving energy 
coincides with ©, up to a constant term ca 2 12w z . 

3.5 Connection to (the nonholonomic failure of) Noether theorem. In the set¬ 
ting of Theorem [5] it is natural to view the conservation of the moving energy 

E*l,c v ,m 0 = El,M o ~ Jv,m 0 ! 

where 

J 9 ,M 0 = Jq | Mo • 
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as related to the invariance of the system under the action T. If the energy El,m 0 is 
not conserved, then the conservation of the moving energy Ef c Mq is only possible if 
J V} m 0 is not conserved. Thus, Theorem [2] produces a conserved quantity from the sum of 
two quantities—the energy and a component of the momentum map—that, at variance 
from what would happen if the system was holonomic, are not conserved. In a way, this 
first integral seems to be produced notwithstanding—or perhaps thanks to—the failure 
in nonholonomic mechanics of two cornerstones of Lagrangian mechanics: conservation 
of energy and Noether theorem. 

It has some interest to understand why, at least in the symmetric case considered 
here, the mechanism of Theorem [Q that obviously has no interest in the Lagrangian 
case, is instead of interest in the nonholonomic case. In the absence of nonholonomic 
constraints (Mq = Mq = TQ ), Theorem [T] is obviously true: the time-independence of 
the two Lagrangians L and L implies that both functions El and Ef c are first integrals 
of the Lagrangian system described by the Lagrangian L. However, since in Lagrangian 
mechanics the momentum map of a lifted action that leaves the Lagrangian invariant is 
conserved, this mechanism can hardly be seen as disclosing a ‘new’ first integral Ef c : 
the difference Ef c — El is a component of the momentum map, and it is thus a first 
integral for the very same reason of symmetry that underpins the possibility of passing 
to moving coordinates without introducing a time-dependence in the Lagrangian. 

Explaining why things are different in the nonholonomic case—and how they are 
different—requires exploiting the role of the reaction forces, along the lines of ilT, 12] . 
The hypothesis of ideality of the constraints assumes that the constraint can—a priori— 
exert all reaction forces that he in (X>o)°. However, the set of reaction forces that is 
actually exerted when the system Xl,q,m 0 is i n a configuration q £ Q with any possible 
velocity q £ (Mo) q is given by 

TZq := [J R L ,M(q,q) 

q£(M 0 )q 

and can be (and typically is) smaller than (Vo)° q . Therefore, the annihilators of 7 Z q can be 
(and typically are) larger than the fibers of ( V t ) q . These annihilators are the fibers of a 
distribution 1Z° on Q , which was introduced in m (in the case of time-independent linear 
constraints, but the generalization to the case of time-independent affine constraints is 
straightforward m and was called the reaction-annihilator distribution. We refer to 
these works for further details and we limit ourselves to note that 

(' Vo) q CK° q \/q,t. 

Proposition 4. [12] Consider a time-independent nonholonomic system with affine con¬ 
straints Xl,q,m 0 > an d 1C Mq = Vq + fo ■ 
i. The energy El,m 0 is conserved if and only iff, o is a section of 1Z°. 

ii. Assume that the Lagrangian L : TQ —> R is invariant under the tangent lift of an 
action of a Lie group G on Q, namely L o 'S?J G = L for all g £ G. Then, for any 
r] £ g, J V: m 0 is a first integral of Xl,q,m 0 if an d on ly if^'v is a section oflZ°. 

Assume, thus, that the energy of a time-independent nonholonomic system Xl,q,m 0 
with affine constraints is not conserved. By Propositions] this happens if and only if £o is 
not a section of 7 Z°. By Theorem^ under Hypotheses (HI) and (H2), the existence of a 
conserved moving energy E’f c Mq depends on the existence of an infinitesimal generator 
Y v such that the difference £o — Y v is a section of Vq. Since T>q C 7 Z° and £o is not a 
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Figure 1: The sphere in a rotating convex surface of revolution. 


section of 1Z °, this necessarily requires that Y v has a nonzero component off 1Z° and, still 
by Proposition |4l the momentum J v ,m 0 is not conserved. Thus, at the basis of the non¬ 
conservation of the energy and of the component of the momentum map—that makes it 
possible for them to add up to give a conserved function—there is the same reason: £o 
is not a section of 1Z°. 

Remark. One might take as well a different point of view, and see the function E £ as 
the momentum map of the action, in the extended phase space, given by a combination 
of time-translation and of the lift of a one-parameter subgroup of T. 

4. Integrability of a sphere rolling on a rotating surface of revo¬ 
lution 

We outline now an application of the existence of a conserved moving energy to the 
class of systems, considered by Borisov, Mamaev and Kilin in [B], that are formed by 
a heavy homogeneous solid sphere constrained to rotate without sliding on a moving 
surface of revolution; specifically, the surface rotates—relatively to an inertial frame— 
with uniform angular velocity $7 around its figure axis, which is assumed to be vertical 
(that is, directed like gravity). Describing the system in an inertial frame {O; x,y, zj 
and using time-independent coordinates, as done in [(Jj, leads to a time-independent 
Lagrangian; we will assume that the 2 -axis coincides with the figure axis of the surface, 
see Figure 1. The case $7 = 0 is classical. Its study goes back to Routh m > who also 
considered special cases with SI ^ 0 (see also [22). When $7 ^ 0 and the surface is a 
horizontal plane, the system reduces to the sphere on the turntable of section 3.1. 

In all cases, either with $7 = 0 or with SI ^ 0, the (time-independent) constraint 
manifold Mq has dimension 8 and the system has an obvious S 1 x SO(3) symmetry. 
Reduction leads to a four-dimensional system which has two independent first integrals, 
that we will denote J\ and J%. We will use the same symbols J\ and J 2 also for the lifts of 
these functions to the unreduced phase space, which are first integrals of the unreduced 
system. The existence of the two first integrals J\ and J 2 , when $7 = 0, was known 
already to Routh [24[ [22]. Their existence when $7 ^ 0 has been proven in [6] • 

The integrability of the case $7 = 0 has been studied by Hermans [16: and Zenkov J7]. 
When $7 = 0 the constraint is linear and the energy £l,m 0 is conserved. As a result, the 
unreduced system has the three first integrals J±, J 2 and £x,m 0 , which are independent in 
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an open subset M { * of the phase space M 0 (specifically, they are independent everywhere 
except on motions in which the center of the sphere either moves horizontally or is at 
rest). Since the function El,m 0 is S 1 x SO(3)-invariant, the reduced system has three 
independent integrals as well. Moreover, if the surface is upward convex and the sphere 
rolls inside it, as in Figure 1, then the common level sets of these three integrals in the 
four-dimensional reduced phase space are compact, and hence are closed curves, and 
the reduced dynamics is periodic [161127] (reference [TB] uses a different argument; for 
details on the properties of the first integrals see BSD- Since the symmetry group is 
compact and acts freely on Mq, this in turn implies that the unreduced dynamics in 
Mq is quasi-periodic on tori of dimension up to three. This was proven in [16] using 
a reconstruction result from periodic dynamics, originally due to Field and Krupa (see 
particularly m m m ® m nsi). Reference m reaches the same conclusion, but 
restricted to the motion of the center of mass, that undergoes quasi-periodic motions on 
tori of dimension up to two. 

When 0 the constraint is affine, not linear, and even if the constraint and the 
Lagrangian are time-independent the energy is not conserved. Therefore, the argument 
used for the case = 0 is not directly applicable. A different approach has been taken 
by Borisov, Mamaev and Kilin, who proved that the reduced system has an invariant 
measure and, using Jacobi theorem [2], deduced from this and from the existence of the 
two first integrals J\ and J 2 that the reduced system is integrable by quadratures [Bj. 

But as we now prove, a conserved ‘moving energy’ exists in this problem, and much 
stronger results can be obtained. Leaving for a future work a detailed study of the 
problem, we limit here ourselves to some conclusions that can be obtained combining 
Theorem |T] with some general arguments (essentially, continuity from the case Q = 0): 

Theorem 3. Consider the system formed by a heavy homogeneous solid sphere that rolls 
without sliding on a surface of revolution, which rotates with constant angular velocity Q 
around its figure axis, aligned with gravity. Then, at least for SI not too large: 

1. The reduced system has three first integrals, which are independent in some open 
nonempty subset of the four-dimensional reduced phase space. 

If, moreover, the surface is upward convex, and the sphere rolls inside it, then: 

2. There is a nonempty open subset of the reduced phase space where the reduced dy¬ 
namics is periodic. 

3. There is a nonempty open subset of the phase space of the unreduced system in which 
motions are quasi-periodic, on tori of dimension up to three. 

Proof. Let r G R 3 be the vector of the coordinates of the center of the sphere relative to 
the considered inertial frame {O; x, y, z }, R G SO(3) be the matrix that fixes the attitude 
of the sphere relatively to that frame, and w 6 i 3 be the angular velocity in space of the 
sphere relative to that frame. 

The holonomic constraint that the sphere is in contact with the surface of revolution 
can be modelled by imposing that the vector r is constrained to a (fixed, time inde¬ 
pendent) surface of revolution E, that we embed in R 3 9 r. The configuration space of 
the holonomic system is thus Q = E x SO(3) 9 (r,R) and the phase space TQ can be 
identified with E x R 2 x SO(3) x R 3 9 (r, r, R,uj). The Lagrangian is the restriction to 
E x M 2 x SO(3) x R 3 of the function L : R 3 x R 3 x SO(3) x R 3 —> R given by 

L(r,r,R,u) = i||r|| 2 + ica 2 ||w|| 2 - gr 3 , 
where the constants have obvious meanings, see ©, and is time-independent. 
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The constraint of rolling without sliding leads to a time-independent nonholonomic 
system with affine constraints, with constraint submanifold an 8-dimensional affine sub¬ 
bundle Mo of S x l 2 x SO (3) x R 3 and Lagrangian the restriction of L to Mq. The 
Lagrangian and the constraint manifold are invariant under the tangent lift of the action 
of S 1 x SO(3) on Q given by 


* {e , P) (r,R) = (H e r,H e RP) (12) 

where Hg is the 3x3 matrix of rotation of angle 6 around the third axis. Once the 
kinematical states with the sphere sitting at the point r = 0 and spinning about the , 2 - 
axis have been removed from phase space, to prevent the need for singular reduction, the 
(regularly) reduced phase space is a 4-dimensional submanifold of E x R 2 x R 3 E3 (r, r, oj). 

We pass now to time-dependent coordinates (s, S ) in Q with 

Q(s,S,t) = (Ha t s,HntS) (13) 

and lift them to a time-dependent coordinate change DC : (s, s, S , u, t) M- (r, r, R , ui, t) in 
TQ. In this coordinate system the surface is at rest; therefore, the constraint of rolling 
without sliding is linear and time-independent and defines a linear subbundle Mq of 
£ x R 2 x SO(3) x R 3 . Due to the symmetry of the system, the Lagrangian L = L o DC 
is time-independent as well. (Incidentally, since Q is the restriction of the action ik to a 
one-parameter subgroup, the time-independence of L follows as well from the argument 
used in the proof of Theorem [2] we mention also that 

L{s,s,S,v) = + Cle z x s || 2 + ^ca 2 ||u + Oe-1 | 2 - gs 3 ). 

Thus, X ^ q m 0 is a time-independent nonholonomic system with linear constraints. More¬ 
over, by item i. of Proposition [3] the function E p c is time-independent. By Theorem [Tj 
we conclude that the original system X[^q m 0 has the conserved moving energy E^ c M . 

We may now prove the three statements. 

1. Since the constraint manifolds Mq and Mo, the Lagrangians L and L and the 
change of coordinates are S' 1 x SO(3)-invariant, the function E^ c Mg has this very same 
invariance property and descends to a first integral of the reduced system. The reduced 
system has therefore the three first integrals J\ , J 2 and E p c Mq . From the expressions 
for J\ and J -2 given in [5] and from the expression above of L, one sees that these integrals 
depend continuously on f1 (for Cl = 0, E £ c Mg reduces to the energy El } m 0 ), and the 
same is obviously true for the constraint manifold Mo- It is known that, when Cl = 0, Ji, 
J 2 and E p c Mq are the components of a submersion from an open nonempty set M reg of 
m 0 to r 3 123 mg. Continuity implies that, at least for Cl sufficiently close to zero, the 
map (Ji, J 2 ,E p c Mo ) is a submersion from an open nonempty subset of M reg to R 3 . 

2. When Cl = 0, the level curves of (Ji, J 2 , E P C Mq ) in M reg are compact [27, 10], 
hence bounded. It follows that, for each Cl sufficiently close to zero, there is an open 
nonempty subset of M reg where the submersion (Ji, J 2 , E^ c Mg ) has bounded (hence 
compact, being a submersion) level sets. 

3. According to the mentioned reconstruction results from reduced periodic dynamics 

by Field and Krupa, if the group is compact and acts freely, then each ‘relative periodic 
orbit’ (that is, the group orbit in the phase space that projects over a periodic orbit of 
the reduced system) is fibered by tori of dimension up to 1 + p, where p is the rank of 
the group, on which motions are quasi-periodic. In our case, p = 2. □ 
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These are clearly partial results, that should be completed and extended under several 
aspects. We shortly indicate some of them. 

First of all, the regions of the phase space where the dynamics is quasi-periodic should 
be identified, and it should be understood how they depend on SI and on the shape of 
the surface. 

A complementary question concerns the behaviour of motions that are not quasi- 
periodic, if present. 

There are also interesting questions about the geometry of the (singular) foliation by 
invariant tori. The proof given above shows that each relative periodic orbit is fibered 
by tori of some dimension between 1 and 3, but it does not ensure that this dimension 
is the same across different relative periodic orbits and that the invariant tori are the 
fibers of a fibration of (an open subset of) the phase space. This property is important, 
because it implies the existence of the appropriate number of first integrals that are 
usually associated to integrability (for some results on this point in the case 0 = 0 see 
0 ). ' 

Finally, even if the reconstruction procedure gives generically a fibration by invariant 
tori of dimension 3, it might happen that for certain shapes of the surface there are 
resonance conditions among the frequencies of all the quasi-periodic motions—or equiv¬ 
alently, there exist additional first integrals—and there is a fibration by invariant tori 
of smaller dimension, either 2 (all motions have two frequencies) or 1 (all motions are 
periodic). An instance of this possibility is met in the limiting case of the sphere on the 
turntable: since the SO(3)-reduced system has periodic dynamics, and SO(3) has rank 1, 
the unreduced motions are quasi-periodic on tori of dimension at most 2. 

This analysis (which has not yet been performed completely for the case f2 = 0, either) 
requires manifestly an approach different from the general one used in this section, and 
will be done elsewhere. 

Remark. The integrability result for the reduced system given in [6], based on Jacobi 
theorem, could be strengthen if the common level sets of the two first integrals J\ and 
J 2 were compact. Under such a hypothesis, Jacobi theorem implies that these level sets 
are two-dimensional tori and motions on them are linear after a time reparameterization. 
However, not only the level sets of J\ and J 2 are unlikely to be compact (these two 
functions are linear in some of the coordinates on the reduced phase space), but because of 
the time reparameterization, this result would be much weaker than those in Theorem [3] 
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